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Abstract 

The sigma models on projective superspaces CP JV+M ~ 1 ' iv with topological angle 9 = n mod 2-k 
flow to non-unitary, logarithmic conformal field theories in the low-energy limit. In this paper, 
we determine the exact spectrum of these theories for all open boundary conditions preserving 
the full global symmetry of the model, generalizing recent work on the particular case M = 
[C. Candu et al, JHEP02(2010)015]. In the sigma model setting, these boundary conditions 
are associated with complex line bundles, and are labelled by an integer, related with the 
exact value of 0. Our approach relies on a spin chain regularization, where the boundary 
conditions now correspond to the introduction of additional edge states. The exact values of 
the exponents then follow from a lengthy algebraic analysis, a reformulation of the spin chain 
in terms of crossing and non-crossing loops (represented as a certain subalgebra of the Brauer 
algebra), and earlier results on the so-called one- and two-boundary Temperley Lieb algebras 
(also known as blob algebras). A remarkable result is that the exponents, in general, turn out 
to be irrational. The case M = 1 has direct applications to the spin quantum Hall effect, which 
will be discussed in a sequel. 
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1 Introduction 

Two-dimensional sigma models on supergroups and supercosets play a fundamental role in several 
areas of theoretical physics, such as phase transitions in 2 + 1 non interacting disordered electronic 
systems [1, 2, 3] (for reviews see [4, 5]), or the AdS/CFT duality [6, 7, 8, 9]. The study of these 
models is on the other hand quite difficult, for a variety of reasons all related to the non unitarity of 
the target. Progress has thus been rather slow, but has picked up pace recently, thanks to the use 
of mini superspace technology [10], a better understanding of algebraic aspects [11, 12, 13, 14], and 
the introduction of lattice regularizations [15, 16, 17, 18, 19, 20, 21] inspired from network models 
[22, 23, 24]. 

A particularly interesting question for critical super sigma models concerns the classification of 
their conformal boundary conditions. While for Wess Zumino Wittcn models the question is well 
understood [25, 26], the situation for other models — such as the superprojective sigma models at 
9 = 7T — is much less under control. In a recent paper [16], the sigma models on 

Cp iV-i|JV = u(JV|JV)/(U(l) x U(JV- l\N)) (1) 

at 9 — 7r were considered, and the set of conformal boundary conditions invariant under the full 
group PU(iV|7V) = XJ(N\N)/XJ(1) classified. These boundary conditions are parametrized by an in- 
teger, and correspond, in string theory parlance, to volume filling branes equipped with a monopole 
line bundle and a connection. In sigma model language, the existence of these boundary conditions 
has to do with the fact that, while the bulk theory only depends on the value of 9 modulo 2ir, when 
there is a boundary, the exact value of the topological angle matters, and enters, for instance, the 
classical equations of motion [27, 28, 29]. 

The sigma model considered in [16] is peculiar because it is critical for a large domain of values of 
the coupling constant g\ including the minisuperspace limit g 2 G — > 0. The solution using harmonic 
analysis on the target in this limit is a precious help in finding out the exact spectrum for all 
values of g 2 . up to the critical value g\ = 1 (in normalizations of [16]) beyond which the sigma 
model exhibits different behavior. For other types of superprojective sigma models, such a line of 
fixed points is not available, and one must directly solve the model at finite coupling. The lattice 
regularizations of [17] constitute a powerful means of obtaining such a solution, and we shall take 
this route in the present paper. 

The main objective of this work is to extend the CP^ -1 '^ study of [16] to the considerably more 
difficult and richer case of CF N+M -^ N . The use of lattice regularizations in this context is not 
only a practical tool, but also provides a physical intuition for understanding conformal boundary 
conditions in terms of edge states. 1 This in turn has applications to the description of the transition 

1 We note that the role of edge states in sigma models has recently been revisited in [30], although the relation 
with our and other's work is not clear to us. 
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between plateaus in the spin quantum Hall effect [31], which corresponds to the M — 1 case of the 
general construction. We will discuss the implications of this observation in a subsequent work [32] . 

Our paper is organized as follows. In section 2, we revisit some of the main results in [16] 
and discuss boundary conditions in sigma models. We review in detail the easy case of M = 0, 
N = 1 (symplectic fermions) which will serve as a benchmark for the subsequent developments. 
In section 3, we discuss briefly the role of boundaries in mapping spin chains to sigma models, 
together with the concept of edge states. In section 4, we discuss the spin chains relevant for the 
solution of the boundary superprojective sigma models, together with their geometrical formulation 
in terms of loops and in terms of lattice algebras. Section 5 is the most technical: this is where we 
obtain the solution of the loop model based on earlier results on one and two boundary Temperley- 
Lieb algebra, together with exact solutions for CP^" 1 '^ and exact diagonalizations of finite size 
systems otherwise. In section 6, we discuss in details the relation between the spectrum of the spin 
chains and loop model, while section 7 contains conclusions. Some particularly technical results are 
discussed in the appendices. 

1.1 Notations 

For the reader's convenience we collect here some notations to be used in this paper. 

• N and N + M denote respectively the number of bosonic and fermionic coordinates of the 
superspaces used. 

• L is the number of pairs of alternating representations in the bulk of the spin chain. 

• m is the number of additional representations on the left boundary of the spin chain. 

• n is the number of additional representations on the right boundary of the spin chain. 

• m = mm(m, n) is the minimal number of uncontractible pairs of boundary representations. 

• C and 1Z are integer (monopolc) numbers labelling different boundary conditions at the ends 
of the world-sheet in the sigma model formulation. 

• (3 is the fugacity of bulk loops in the geometrical formulation of the model. 

• P\ and @2 are the fugacities of loops touching respectively the left and the right boundary in 
the two-boundary loop model, and j3\2 is the fugacity of loops touching both. 

• £>2L is the Brauer algebra on 2L strands. 

• ■AiL,m,n is the subalgebra of the Brauer algebra involved in the formulation of the lattice 
model. 
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• j is the number of non contracted pairs of bulk lines. 

• k is the number of non contracted pairs of boundary lines. 



• p is the minimal model index, corresponding to setting (5 — 2 cos 




• h rtS is the Kac table of conformal weights. 

2 Boundary sigma models 

We first briefly review the definition of the £^> N + M - 1 \ N mo dels, and then move on to describe the 
role of boundaries in the sigma models. 

2.1 Superprojective sigma models at 9 = n 

Complex projective superspaces (£]p N + M ~ 1 \ N are built much like their bosonic cousins [17, 16]. 
Begin with superspace C N +M\N. the N + M complex bosonic coordinates are denoted by z a and 
we use £ a for the N fermionic directions. Within this complex superspace, consider the odd (real) 
dimensional supersphere defined by the equation 

N+M N 
a—l a—1 

The supersphere S 2N+2M ~ 1 \ 2N carries an action of U(l) by simultaneous phase rotations of all 
bosonic and fermionic coordinates, 

Za * ^ Z a • £a * ^ £a ■ (3) 

Note that this transformation indeed leaves the constraint invariant. The complex projective su- 
perspace £P N+M - m is the quotient space ^ 2JV+2M - 1 l 2Ar /C/(l). 

Functions on the supersphere 5' 2Ar + 2M - 1 l 2Ar carry an action of the Lie supergroup \J(N + M \N). 
These transformations include the phase rotations (3) which act trivially on £ff N + M ^ 1 \ N _ Hence, 
the stabilizer subalgebra of a point on the projective superspace is given by u(l) x u(N + M — 1\N) 
where the first factor corresponds to the action (3). We conclude that 

cf n+m-i\n = u(JV + Af |JV)/(U(1) x XJ(N + M- 1\N)) . (4) 

Their simplest representative of interest to us is CP ' 1 i.e. the space with just two real fermionic 
coordinates. The sigma model with this target space is equivalent to the theory of two symplectic 
fcrmions, which has been extensively investigated, as for example in [33, 34]. 

The construction of the sigma model on £^> N + M - 1 \ N closely parallels this geometric con- 
struction. The model involves a field multiplet Z a = Z a (z,z) with N + M bosonic components 
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Z a = z a ,a = 1, .. . ,N + M and TV fermionic fields Z a = £ a - N -M,a = N + M+ 1, . . . , 2N + M + 1. 
To distinguish between bosons and fermions we introduce from now on a grading function | • | , which 
is when evaluated on the labels of bosonic and 1 on the labels of fermionic quantities. In addition 
we also need a non-dynamical U(l) gauge held a. With this held content, the action takes the form 
(with a summation convention on the index a) 

s = 27 / d2z ( 9 n - m ^ z i^ + m ^) z - ~^J d2z ^ d » a » ( 5 ) 

and the helds Z a are subject to the constraint Z\Z a = 1? The integration over the abelian gauge 
held can be performed explicitly and it leads to the replacement 

«m = \ [Z^Z a - (d,Zl)Z a ] . (6) 

The term multiplied by 9 does not contribute to the equations of motion for a^. As its bosonic 
counterpart, the Cp Ar + J1/_1 l Ar sigma model on a closed surface possesses instanton solutions. The 
corresponding instanton number is computed by the term — which we will denote by Q and refer to 
as the 'topological term' — that multiplies the parameter 9. Since Q is integer- valued, the parameter 
9 = 9 + 2tt can be considered periodic as long as the world-sheet has no boundary. 

The target supermanifold being a symmetric superspace, the metric on the target space is unique 
up to a constant factor, so g 2 a and 9 are the only coupling constants. The perturbative beta function 
is the same as the one for CP M_1 [35] 

^- = P{gl) = Mgi + 0{gl) (7) 

The beta function for 9 is zero in perturbation theory, and that for g^ is independent of 9. 

For M > 0, the coupling is weak at short length scales, but flows to strong values at large 
length scales. For 9 ^ it (mod 2w) the U(iV + M \N) symmetry is eventually restored, and the 
theory is massive. When 9 = it (mod 2ir) and M < 2, the model flows to a non-trivial flxed point. 
The corresponding bulk conformal held theory has been studied in [17, 19], and presents many 
fascinating features. The main purpose of this paper is to study its boundary properties in more 
details. 

For M = 2 and N = we recover the usual flow in the 0(3) sigma model; the conformal field 
theory in that case is the SU(2) level 1 Wess Zumino model. For M — 0, the beta function is in 
fact exactly zero, and the sigma model exhibits a line of fixed points [17, 16, 36, 4, 37]. The case 
M = 1 is particularly interesting, because it is related with percolation, and with the spin quantum 
Hall effect [31, 32]. 

2 Note that we eliminated the radius p of the complex projective space in favor of a coupling g^ 2 entering the 
action in front of the metric. Equivalently, we can set g 2 = 1 and work with a radius parameter p appearing in the 
modified constraint ZaZ a = 4p 2 . 
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2.2 The role of the boundaries 

While the model on a compact oriented manifold exhibits properties that do not depend on the 
exact value of 9 provided 9 — n (mod 2tt) , things are bound to be quite different in the presence of 
a boundary, as first noticed in [29]. Strictly speaking, in a finite system with boundaries, there are 
no well defined topologically distinct sectors. One does not expect the behavior of the RG for the 
bulk properties to change — the standard analysis should still apply provided the equivalent of the 
instantons (configurations such that the 'topological term' Q approaches a constant in the core, far 
from the boundary) are well localized within the system. But boundary properties are expected to 
now depend on the exact value of 9; in particular, 9 can now be expected to affect perturbation 
theory. 

To make things more concrete, we restrict to sigma models on the strip £ = [0, tt] x K or, 
cquivalcntly by conformal transformation, the upper half plane z — x + iy, y > 0. The boundary 
conditions induced by the 'topological term' are 

Q 

(d y + ia y )Z a = -gl(d x + ia x )Z a , 

(d y - ia y )Zl = g 2 a (d x - ta x )Zl 

for z — z < and a similar condition along the right half z = z > of the boundary. While 
in the weak coupling limit the value of 9 is irrelevant (and we recover purely Neumann boundary 
conditions), in general, 9 appears explicitly. 

Of course, except in the special case M = [17], we will need 9 = 7r(mod 2tt) to have a conformal 
field theory in the low energy limit, but we can now expect properties to depend on the exact value 
of 9 itself. 

It was indeed argued in [16] that a more general family of \J(N + M \N) symmetric boundary 
conditions can be obtained which are expressed through the conditions 

(d y + ia y )Z a = Q^Kdx + ia x )Z a , 

(9) 

(d y -ia y )Zl= -B ig l(d x - ia x )Zl 

for z — z < and a similar condition with Oi replaced by O2, along the right half z = z > of 
the boundary. The parameters 81 = 2£ + 9/n and 62 = 21Z + 9/ir involve now integer (monopole) 
numbers £,1Z. 

In the case M = 0, it was shown in [16] how, to every choice of integer 1Z = C, is associated a 
conformally invariant boundary condition. The case 1Z =/= C then corresponds to the insertion of a 
boundary condition changing operator. 

For < M < 2, we expect similarly that conformally invariant boundary conditions at the 
fixed point theory are induced by these boundary terms. In other words, we expect that, in the 
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conformally invariant fixed points of superprojective sigma models C¥ N+M - llN at 9 = it, there is 
a discrete family of conformal boundary conditions labelled by the integer R £ Z: these are the 
subject of our study. 

We note that the only dependence on the exact value of 9 comes from the boundary conditions. 
A model with a given value of 9, £, 1Z is identical with the model with 9 + 2pn 7 C — p, 1Z — p, for any 

pel,. 

2.3 The case CP D|1 

To make things more concrete, we note that the simplest case of CP ' 1 (N = 1, M = 0) is equivalent 
to symplectic fermions. With the bulk action 

S=^Jd 2 zd^ (10) 

the boundary conditions are then of the form 

d„? = -Qg*d x ? (11) 

They can be represented in terms of a glueing automorphism 

$)-»($) ^ 

where 

If we now have two different values of on the left and right boundaries, going around the insertion 
of the boundary condition changing operator in the complex plane gives rise to a monodromy 
expressed by 

_ pp -i_ n + w(e 1 -e 2 ) 
° 12 - fil ° 2 -«-w(e 1 -e 2 ) (14) 

where we have set k = 1 + ,g*0i0 2 . Of course, fii 2 is of the form 

Oi2 = [ o e -L J (15) 
with the twist parameter A given by 2cos2-7rA = Tr f2 12 : 

rnM (i + g^e 1 e 2 ) 2 -(9 1 -e 2 ) 2 g^ 

The twist parameter vanishes when C = 1Z. In general, the ground state of the theory (which has 
central charge c = —2) scales with the conformal weight 

fcf = ^A(A-l) (17) 
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and the full operator content in this sector is encoded in the U(l) x Vir character [33] 



dfj,,x = tr ( e 



,2ifj,.Jo Lo-c/24 



) 



2i7TLlX 




,n+A-l 



)(l + c 



— 2i7Tfl 7l — X 



) 



(18) 



71=1 



or 



1 -2i7T 



juA ^ ' e 2i7T^im^i(m+A-i) 2 



(19) 



where 9 denotes the modular parameter and rj is the Dedekind eta function. 

The problem is thus fully solved in this case. In this paper, we want to address the considerably 
more difficult case of M ^ 0, in particular the case M = 1 which is relevant to the spin quantum 
Hall effect [31, 32]. 

3 Boundary conditions in super spin chains 
3.1 Super spin chains 

As mentioned in the introduction, our strategy to solve the boundary sigma model is to use a lattice 
regularization. There is a well known, profound relationship between sigma models and spin chains, 
going back to the earliest developments in the 0(3) case [38]. In fact, the conformal fixed points 
in the bulk £jp N + M - 1 \ N sigma models at 9 = tt have been studied in [17, 19] using homogeneous 
super spin chains that represent the strong coupling region. The mapping to the sigma models 
follows the well known argument in non graded [39, 40] as well as graded [41, 31] cases. 

The simplest chain (which will be described in considerably more details below) is obtained by 
alternating the fundamental representation V of sl(N + M\N) and its conjugate V* , and choosing 
an antiferromagnetic Heisenberg coupling between nearest neighbours (since the product V ® V* 
decomposes generically on two representations — the identity and the adjoint — this is the most 
general nearest neighbour interaction). This coupling in turn can be conveniently recognized as a 
representation E sl of the Temperley-Lieb algebra acting on (V <E) V*)® L , and the Hamiltonian with 
open (free) boundary conditions can be written as 



i=l 

Detailed expressions in terms of the natural vector basis of V, V* are given in [17, 18] and we do 
not reproduce them for now (but see below). An important point for further study is that the 



2L-1 




(20) 
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generators E? satisfy the Temper ley-Lieb relations [42] 

(Eff = ME? 

[E?,Ef] = 0, for |i - J] > 2 . (21) 

Instead of the Hamiltonian, it is sometimes useful to use a transfer matrix, the evolution operator of 
the ID chain in discrete imaginary time. The choice of alternating representations in the spin chain 
corresponds then to having oriented edges carrying respectively V and V* as illustrated in figure 1 , 
making the potential relationship with network models of quantum localisation more transparent. 
The transfer matrix itself is then built as a product of two diagonal-to-diagonal transfer matrices 
associated to two consecutive layers of the lattice, which evolves the states one unit in time. 




Figure 1: The two dimensional lattice corresponding to the alternating super spin chain. Edges 
carrying the fundamental representation have up arrows while those with the dual down arrows. 
The transfer matrix is a product of two transfer matrices acting on two consecutive layers T = XY , 
imaginary time flows from bottom to top. 

For < M < 2, the model thus defined flows to the fixed points discussed in [19]. The 
Hamiltonian with periodic boundary conditions has been studied in [17]. When M = 0, recall 
[17, 16, 36, 4, 37] that the sigma model admits a line of fixed points parametrized by g 2 : the model 
we have just defined corresponds to a particular value of g 2 a (= 1). Other values can be obtained 
by adding an exactly marginal interaction corresponding to the exchange of next nearest neighbors 
degrees of freedom V or V* (since the products V <g> V or V* ® V* decompose generically on two 
representations, this is the most generic nearest neighbor coupling). In the Hamiltonian language, 
this means 

2L-1 2L-2 

H = ~ E E ?~ w E ( 22 ) 
z=l i=l 

where -Pj S ' +2 is the operator permuting states in sites i and i + 2. Note that the Pf] +2 t erms do not 
break the sl(N + M\N) symmetry [36], i.e., it does not mix up the V and V* spaces. A proposal 
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for obtaining the conformal boundary conditions starting from this spin chain was then made in 
[16]. There, a family of four possible chains was considered, corresponding to the following spaces 
(with factors labelled from to 2L + m + n — 1) and Hamiltonians 

V® m <g> (V <g> V*)® L <g> : H vv * = fl£ ft + H hulk + HY* ht 



V® m ® (1/ ® V*)® L ® 1/®" : F yv = #£ ft + i?bulk + ffj 



'loft "r ■'■'bulk T ■'■'right 

(y*j(8m g (V ® T/*)® L ® V"®" : H v ' v = + i? bulk + F r Y ght , 



(23) 



where the bulk Hamiltonian is 



2L+m-2 2!+m-3 

ffbuik = - E £f - w £ i# +a > ( 24 ) 



7.=m 



while the boundary Hamiltonians are as follows 

m—X 2L+m+ro— 2 

<ft = - u E P ti+i bright = -« E P M+i (25) 

1=0 i=2L+ro-l 
m-2 

= - « E P W - W ' P *Um+l - t'Kll (26) 
i=0 

2L+m+n-2 

bright = — ^"-^2L+m-l — w " P2L+m-2,2L+m — u E ^M+l ' (^"0 

i— 2L+m 

In words, the Temperley-Lieb generators operate only in the bulk spaces, whereas the permutation 
operators act on the boundary spaces and on the first physical space of the same type as the 
boundary spaces. 

Numerical evidence was presented [16] to the effect that these chains correspond to the sigma 
model with monopolc boundary conditions according to the correspondence 

V® m ® {V®V*)® L ® (V*)® n : C = +m TZ = +n , (28) 

V® m ® {V®V*)® L ® V® n : C = +m TZ = -n, (29) 



c = 


+m 


TZ = 


+n 


c = 


+m 


TZ = 


—n 


c = 


— m 


TZ = 


+n 


c = 


—m 


TZ = 


—n 



(V*'f m <g> [V ® V*)® L ® (V*)® n : C = -m TZ = +n, (30) 

(V*f m ®(V®V*)® L (g> V® n : C = -m TZ = -n . (31) 
It was also argued that 

6 = 7T (32) 
and that, for the value w = to which we restrict here, 

M = 0, g 2 a = l. (33) 
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3.2 Mapping on sigma models: the role of boundaries 

We can give a quick heuristic derivation of these results following [43, 44]. Consider the standard 
mapping of the spin chain onto the sigma model. Calling the elementary Berry phase O(i), the 
total phase for the usual antiferromagnetic chain of spin s would read 

$ to P = s V(-i)<n(i) « - / — dx = - pj(2L) - n(o) + 4^g] (34) 

where Q is an integer. Let us now suppose we have for instance the first chain in equation (28). 
Since on the boundary we have representations of the same kind, the Berry phase on that side 
comes unstaggered, so we have (from now on we put s = \, which is the right value for our spin 
chain [17]) 

2L 

$ top = -|fi(0) + £fi(2L) + \ * ^^(2L) - + ttQ (35) 

i=l 

giving immediately rise to the equations (9) with 9 — ir. Meanwhile, say for the second chain in 
(28) we have 

$ top = -|n (Q ) _ 2n(2i) + 1 ^(-l)^(i) « ~ n + 1/2 ft(2L) - HL+VH n (0) + ttQ (36) 

i=l 

again in agreement with equations (9) this time with C = m but 1Z = —n. Note that in this kind 
of argument, the exact nature of the couplings on the boundary is irrelevant. Note also that it 
does not make a difference whether the spins on the boundary are actually projected onto the fully 
symmetric representation or not: we will discuss this in more details below. 

We thus see that adjusting the £,71 terms in the formal boundary action of the sigma model 
is equivalent to adding extra spins — or edge states — to the boundaries of the spin chain. As com- 
mented earlier, the model with 9 — ir and £ = 7l=P7^0is equivalent to the same model with 
9 = 7r + 2irp and C = TZ = 0. Hence shifting 9 by multiples of 2tt leads to the apparition of 
extra edge states. This phenomenon is well known in the language of QED (see [45]), to which 
the CP M_1 model is equivalent at large M, with M flavors and a weak gauge coupling e 2 ~ jj, 
as reviewed in [39]. The topological term is equivalent to a background electrostatic field F in 
the one-dimensional universe, with F = If this field is too large, it is energetically favorable 
to produce quark-antiquark pairs. For one such pair, I being the distance between the quark and 
antiquark, the difference in energy between the state with and without the pair is 

AE = I [{F±ef - F 2 } 

While it is not energetically favourable for the vacuum to produce a pair if |F| < |e, it becomes so 
if \F\ > |e. Pairs will in fact be produced until F is brought down to a value [i^creenedl < So 
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if 7r < 9 < 37r, one pair is produced, and more generally if (2n — 1)tt < 9 < (2n + 1)tt, n pairs are 
produced. 

This is in fact for F positive. If F is negative, that is if 9 < 0, things are quite similar, only 
one produces this time antiquark-quark pairs. The picture is thus in agreement with the previous 
findings in the case where there are extra representations V p on the left and (V*) p on the right 
(£ = K=p), or {V*) p on the left and V p on the right (C = TZ = -p), with 9 = § + 2TZ = 1 + 2TZ. 
This is once we have identified the chain with no extra representations as corresponding to 9 = n 
(,s=i). 

Restricting now to the case C — 1Z, we observe that another way to have the same physics would 
be, instead of taking a spin chain with s = | in the bulk and adding edge states by hand, to take 
directly a chain with spin s = \ + C = 1+ 2 2£ . In other words, with a spin s (s half an odd integer) 
chain in the bulk, we expect physics of edge states with C = s — h , which corresponds in turn to a 
spin on the boundary with value s' = s ~^ 2 ■ 

In the usual case of the CP 1 model — that is, the 0(3) model and the XXX spin chain — all 
this discussion is somewhat irrelevant. Although edge states still present interesting features, we 
anticipate that adding extra spins (maybe in a higher dimensional representation) on the boundary 
does not give rise to new exponents, and the spectrum will in the end coincide with that of the XXX 
chain with an odd or an even number of sites. In the case of cp Ar + M_1 l Ar sigma models however, 
an infinity of conformal boundary conditions is available, leading to a much richer behavior. 

4 Spin chains and loop models 

In what follows we will derive results for super spin chains of the type described in section 3.1 when 
M is not fixed to zero. In particular, we focus now on the case of the chain in equation (28), which 
we recall: 

V® m ® (V ® V*f L ® (V*)®" , (37) 

where as before V and V* are two graded vector spaces of even dimension N + M and odd dimension 
N, taken respectively as the fundamental and the dual representation of sl(N + M\N). 

4.1 Hamiltonians and transfer matrices 

We give first the explicit expression of the action of operators on this spin chain. Let ej, with 
j = 1, . . . , 2N + M, be a basis of V, and e J , with j = 1, . . . , 2N + M, be the corresponding dual 
basis of V*, and let |j| denote the grade of ej. Then the graded permutation of index 1, with 
< i < m — 1 or to + 2L — 1 < i < m + 2L + n — 2, acts on two copies of the same representation 
by interchanging them with a minus sign if both the elements are odd: 

Pi,l+i e io ® ■ • • ® e jm+2L+n - 1 = (-) |j '* l|j ' < + l| e jo ® • • • <g) e ji+t <g> e k ® • • • ® e J " l+2I '+"- 1 , (38) 
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The Temperley-Lieb operator of index i, m < i < m + 2L — 2 projects onto the singlet in the 
decomposition of V ® V* (and ® V): 



JV+Af 



- 1 = (-)^ii^i^ 1 £ 



fc=l 



e fe <X < 



The model we consider is given by the following Hamiltonian: 

m— 1 2Z+m-2 2L+m+n-2 



.Jm + 2i + n-l 



(39) 



(40) 



i=2L+m-l 



One can naturally associate to this Hamiltonian a vertex model defined on a strip of length 
m + 2L + n, see figure 2. Like before, the orientations of edges are fixed by the presence of V and 
V*, so there are four types of vertices. 




m 



2L 



Figure 2: Vertex model corresponding to the super spin chain with extra edges at the boundaries. 
Here m = n = 2. 



The transfer matrix at the critical point will be of the form T = XY: 

i~(m) L-l+i~(m) L-l+i~(m)+i + (n) 

x= f[(i+ u p£ +lt2i+2 ) n n (i+^w 2 ) 

i=0 i=i(m) i=L+i _ (m) 

i+(m)-l L-l+i(m) L-l+i(m)+i+(ri) 

f= n (i+^+i) n n a+^+i), 

i=0 i=i+(m) i=L+i(m) 



(41) 



where i(x) and i ± (x) are the integer part of x/2 and (a; ± l)/2 respectively. The parameters u,v 
are not necessarily the same as those in the Hamiltonian, but they conceal the same physics: they 
control the strength of the non-trivial part of the boundary interaction. A more general transfer 
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matrix could be introduced with the Ef l interaction chosen different for even and odd i, as is the 
case for the purely alternating model [17]. 



4.2 Formulation within the Brauer algebra 

As mentioned earlier, the generators E 5> obey the defining relations of the Temperley-Lieb algebra. 
Adding the P 5 ' generators as needed to define the Hamiltonian (40), the Temperley-Lieb relations 
are completed by additional ones, which define a certain subalgebra of the Brauer algebra that we 
will denote A2L,m,n- 

It is useful at this stage to remind the reader that the full Brauer algebra on 2L sites £>2l(/3) is 
generated by the permutation operators Pj(:= P^j+x), with i = 0, . . . , 2L — 2, and the Temperley- 
Lieb generators Ei, with i = 0, . . . , 2L — 2, and is specified by another parameter /3 which we will 
call — for reasons that will soon become obvious — the loop fugacity /3. 

It is well known that words in 2?2l(/3) can be represented graphically using diagrams composed 
by rows of dots connected in pairs. The generators / (identity), Ei (Temperley-Lieb operator) 
and Pi (permutation operator) are shown in figure 3 in this graphical representation. The Brauer 
algebra is then the C-span of diagrams thus obtained. The product of diagrams d\ ■ d,2 is defined 
by placing d\ over di and identifying the bottom dots of d\ with the top dots of di, replacing every 
loop formed by the (complex) weight f3. For an abstract introduction to the Brauer algebra, see 
[46, 47], and for a study in the context of super spin chains, see [15]. 




Ei 



— X 5 


5 


e *— 


— * 3 


i — >C—k — 3 


i *— 



Pi = 




Figure 3: Identity, the Temperley-Lieb generator Ei and the permutation operator Pi, in the 
diagrammatic representation. 



The abstract relations defining BiiXfi) are the following: 



P? = l, E?=/3Ei, E i P i = P i E i 



E-, 



PP±iP = P±iPP±i, EiE i±1 Ei = E u 
PiE i±1 Ei = P i±1 Ei, EiE i±1 Pi = EiP i±1 . 



(42) 
(43) 
(44) 
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and for non-adjacent sites, \i — j\ > 1, all generators commute: 



P P — P P 



EiEj = EjEi, 



(45) 



In the case we wish to study, the full Brauer algebra B 2 l{P) is not needed, since we restrict to 
having permutations on boundary sites only (the Temper ley-Lieb generators still act in the bulk). 
The corresponding subalgebra will be denoted A2L,7n,n((3), and its generators satisfy the relations 



where we note that the two relations (44) have been replaced by the single relation (48) (which 
is implied by the former two within the full Brauer algebra). One can easily verify these defining 
relations by expressing the generators as diagrams. 

We shall not devote much time to studying the algebra A2L,m,n{P) abstractly, as our only 
purpose is to use it as a tool to diagonalize our Hamiltonians. Apart from providing a convenient 
language and intuition, the algebraic point of view indeed will allow us to separate the spectrum 
into different sectors [48], and, within each sector, to considerably reduce the size of the Hilbert 
space necessary for numerical determination or verification of the exponents. 

4.3 Spin chain representation and quotients 

In the spin chain, a representation of the algebra A2L,m,n{P), with loop fugacity /3 equal to STr 1 = 
M, is obtained if we take the action of Pi and Ei given by P/ [ and Ef (see eqs. (38) and (39)). 
This representation is however not faithful: in addition to (46)-(49), the generators acting on the 
spin chain satisfy additional relations. These additional relations define a quotient of A2L,m,n(M), 
which we now describe. 

First we note that the algebra A2L,m,n{M) as defined by (46)-(49) is infinite dimensional. 
We can understand why by considering the simple system with L = m = n = 1. If we call 
W — E1P2PQE1, then the word W p (see figure 4) cannot be reduced using the defining relations, 
and increasing the power p, the dimension of the algebra becomes arbitrarily large. 

It is also clear that this feature is present only if both m and n are non zero, and the larger m 
and n, the larger the number of words that cannot be reduced. A similar situation appears also 
for the two-boundary Temperley-Lieb algebra [49], which indeed is related to our model as will be 
explained later. 

The proper algebraic setting for our model is a quotient of A2L,m,n{M). We will now give the 
expression for the additional relations present in the spin chain representation. Call m = min(m, n) 



Pf = 1, Ef = (3Ei, EiPi = PiEi = E { 



(46) 
(47) 
(48) 
(49) 



PiP i±1 Pi = P i ±iP i P i ±i , EiE i±1 Ei = Ei 



EiPi±\Ei — Ei 



I = I', I',- = EjEi, i:,P, = PjEi , 
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W = ^X^^X^ - W 2 = M ^^ ^^^ ^ W 3 = M 2 



I x X : ^""^x : X X! 



Figure 4: The diagrams corresponding to the words W = EiP 2 PoEi, W 2 and W 3 in -A 2 L,m,n(-^0 
with L = m = n = 1. Here and in the following dotted lines separate bulk and boundary sites. 

and 

e = n Ei n Ei ' ( 5 °) 

where I\ (I2) is the set of even (odd) indices of Temper ley-Lieb generators if m is odd or the set of 
odd (even) indices of Temperley-Lieb generators if m is even. To proceed, let us define the words 
W r by the recurrence relations: 

W = E (51) 

W r — EP m -i ■ ■ ■ P m -rPm+2L-l ' ' ' Pm+2L+r-2W r -\ , (52) 

where r — 1, . . . , m. Also define 

K ir = P m -iWr - P m +2L+i-2~W r , i = l,...,r, (53) 

for r — 1, . . . , m, and call the ideal generated by these elements /. We can now take the quotient 
Q.2L,m,n(M) = A 2 l ,m,n(M) / 1 : we claim that the super spin chain representation of Q 2 L,m,n(M) is 
faithful for a sufficiently large number of states M. 

The mechanism underlying the relations in the quotient can be understood from simple examples 
and easily generalized. Take again the system L = m = n = 1. Then the additional relation in the 
spin chain is P0E1P0P2E1 — P2E1P0P2E1. Indeed If we compute the action of E1P0P2E1 on the 
spin chain, we have: 

M+N 

'./ " io 



E 1 P P 2 E 1 e io ® e h ® ® e J ' 3 = V (_)l*l(liol+li»l) 



(54) 

• efc ® e r ® e r ® e fc , 

that is invariant under the action of Po-f^- 

It is useful to represent diagrammatically the relations. In hgure 5 we have represented them 
for the case L = 3, m = n = 2. 

Note that the additional relations just introduced correctly reduce those words that made the 
algebra infinite dimensional earlier. 

From the operational point of view which will be adopted later for numerical diagonalization, 
the quotient tells us that when intertwined top arcs intersect bottom arcs, their relative intertwining 
does not matter, and could be disentangled. 
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E — E3E5E2E4EQ = 



Wi = EPxP-jE = 



W 2 = EP 1 P Q P 7 P & W 1 




Kn = P1W1 - P 7 W X 



K l2 = PiW 2 - P 7 W 2 = 



K22 = P W 2 - P 8 W 2 = 





Figure 5: Ku = 0,Ki 2 — 0,K 22 — are the additional relations present in the super spin repre- 
sentation for a system with L = 3,m = 2,n = 2. 
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4.4 Geometrical formulation: the loop model representation 

The diagram representation of the generators in our algebra naturally leads to a loop reformulation 
of the model: we can simply expand the product (41) in the vertex model transfer matrix and 
translate each term to obtain a sum over configurations of oriented loops. When viewed along 
the time direction, bulk vertices correspond to the junction of two incoming lines with different 
orientations. These can be contracted (by Ei) or go through (by I), as shown in the lower two lines 
of figure 6. On the other hand, boundary vertices correspond to incoming lines having identical 
orientations. These can cross (by Pi) or go through (by I), as shown in the upper two lines of figure 
6. 

There is thus a total of four possible vertices, which can each be resolved in two ways. 




Figure 6: Vertices of the isotropic vertex model decomposed into the oriented loop configurations. 
u and v are the weights of permutations of respectively up and down links. 

Expanding all the vertices into pieces of oriented loops we end up with a sum over configurations 
like the one illustrated in figure 7. 

Note that while the diagrams used earlier to define the Brauer algebra involve unoriented lines, 
we have here defined the loop configurations in terms of oriented lines. The subalgebra A2L,m,n 
respects this orientation, in the sense that the associated Hamiltonian (40) and transfer matrix 
(41) only generate diagrams in which each loop has a definite orientation. We stress that this 
orientation is a property of the underlying lattice, not of the loop itself. In particular, for a loop 
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Figure 7: A configuration of the supersymmetric lattice model in the case m = n = 2, 2L = 6. 
The loop orientation follows the presence of the fundamental representation of the Lie superalgebra 
sl(N + M\N) (links oriented up) or the dual representation (links oriented down). Each loop carries 
a weight M. 

of a given position and shape, no sum over orientations is implied, and the states entering the 
loop representation do not contain any orientational information. The weight of the loop (3 is a 
parameter of the algebra. 

We shall find it convenient to modify the loop model in such a way that its spectrum is as close 
as possible to the spectrum of the spin chain. To this end we obviously set j3 = M . Moreover, we 
need to go to the quotient Q2L,m,n(M). As a concrete example, the word corresponding to figure 
7 in the algebra A2L,m,n(M), will be reduced as showed in figure 8 using the additional relations 
of the quotient discussed above. 




Figure 8: The two identified words in the quotient Q.2L,m,n(M) corresponding to figure 7. 



These additional rules complete the definition of our loop model, and we can finally relate the 
spectrum of the Hamiltonian (40) in the loop and super spin representations. This exercise is 
similar to the better known cases of the Temperley-Lieb [17] and Brauer [15] algebra. For the 
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Tcmperley-Lieb case, the spectrum of the Hamiltonian (on a given number of sites 2L) acting on 
the sl(N + M\N) chain is the same for every integer TV, the only thing that changes are (non zero) 
multiplicities, that grow very fast with N. In our case, the situation is a bit different. First, we have 
extensive numerical evidence that, calling E(iV) the set of eigenvalues of the sl(N\N) Hamiltonian 
(up to non zero degeneracies), the following relations holds: 

E(l) C E(2) = E(3) =••• = E(iV) = ... (55) 

Moreover, E(2) is also equal to the set of eigenvalues of the Hamiltonian in the loop representation: 
in other words, all spin chains but the st(l|l) one have the same spectrum as the geometrical model. 
The spectrum of st(l 1 1) meanwhile is only a subset; this is discussed in more details in the appendix 
A. 

When we look at sl(N + M\N), M > 0, the numerical evidence is that the spectra are all the 
same (again, up to non zero degeneracies). It is very likely that spin chains for a number of states 
per site greater or equal to 3 provide faithful representations of the quotient algebra Q2L,m,n(M). 
This statement is in line with the findings of [50], which should be possible to generalize in our 
case. 

A final advantage of the loop representation is that we can now study the geometrical model 
for arbitrary values of f3, which is now a simple parameter appearing as the fugacity of loops. 
Note that, except when (3 = M is an integer, we do not know of any spin chain representation of 
the corresponding algebras, be they A, Q or Brauer. This situation has to be contrasted with the 
ordinary Temperley-Lieb algebra, where the XXZ spin chain provides such a representation for all 
values of the loop fugacity. Clearly the presence of loop crossings — even if only between boundary 
lines — makes the problem considerably more involved. These crossings also make impossible a 
Coulomb gas approach to the underlying conformal field theory (for essentially the same technical 
reasons). 

We remark that the bulk part of the Hamiltonian (40) does not follow from the usual Yang- 
Baxter construction for alternating fundamental and dual representations of sl(N + M\N) [51]. 
However, in the bulk, (40) is solvable exactly, since it belongs to the Temperley-Lieb algebra, and 
the realization of this algebra in the I7 g (s[(2))-invariant XXZ spin chain leads to Bethe ansatz 
equations etc (see also [52] for a related situation). We do not know what the situation for our 
model with boundary interactions might be. But it is interesting to observe that in the particular 
case of sl(2), solvable spin chain Hamiltonians with boundary impurity spins can be obtained from 
solutions of the Yang-Baxter equation [53] . More work is needed to clarify the situation. 
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5 Critical exponents 

We now launch into the very technical problem of determining the critical exponents for the geo- 
metrical model. This will allow us later to reconstruct the spectrum of the super spin chains by 
the proper combination of sectors. 

Our loop model in the bulk is nothing but the well known dense loop gas going back to early 
studies of the Potts model and the dense 0(n) model [54, 55]. It is critical for < M < 2. 

The central charge of the dense loop gas is 

C = 1 "K?TT)' < 56 > 

if we parametrize /3 = 2 cos f J^prj ■ We also recall the Kac formula 

_ ((p+l)r-p S f-l 

which is often a convenient means of stating our results on critical exponents. In the following we 
are going to diagonalize the Hamiltonian (40), although as already stressed, we believe that the 
same results can be obtained studying the transfer matrix (41). Conformal invariance relates the 
exponents h to the eigenvalues £i of the Hamiltonian acting on a chain of length L via the finite 
size scaling formula [56, 57]: 

E i = Lf b + 2f s - V ^ + o(J^ . (58) 

In this formula fb and f s are the bulk and surface free energy, c c fr = c— 24ft, and v s = (p+1) sin f^i) 
is the so-called velocity of sound. 

We briefly review how critical exponents are related to the representation theory of the loop 
model (see [58] for more details). On the lattice the geometrical transfer matrix and Hamiltonians 
have a blockwise lower-triangular structure since the Ei annihilate non-contractible lines (or strings) 
in pairs but cannot create any string and Pi does not alter the number of strings. The restriction 
to a sector is done by setting the action of TL generators on two strings equal to zero. Further 
we note that blocks only differing by the configuration of the bottom row of states are identical. 
So for the purpose of computing the spectrum, we can consider reduced states in which only the 
top part of a state on which the Hamiltonian acts is taken. The critical exponents we are going to 
compute, called in the literature watermelon or 2j-leg exponents, are the lowest conformal weights 
in a sector with fixed number of non-contractible lines. 

Before moving on, it is useful to recall some results about the two-boundary loop model (2BLM) 
[49, 59] which we will use later. This model is a dense loop model on a strip in which lines touching 
a boundary get marked in different ways, and marked loops get a weight different from the weight 
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j3 of bulk loops. These boundary weights are fa (resp. fa) for loops touching the left (resp. right) 
boundary only, and fa 2 for loops touching both boundaries. Using the parametrization 



sin((r 1 + l)^) sin (( r2 + l)^_ 

Pi — 7 \ , fa = 



sinfri^-r^ sin (Vo-tt^ 

V P +1 ) \ P +1 ) /.ON 

^ ^ sin [{n + r 2 + 1 - ri 2 )^T) ) sm (fa + r 2 + 1 + ^^TI) ) 

sin ( r w) sin ( r 2iTl) 

it was found that in the sector with no non-contractible lines, the conformal weights appearing are 
[59] 

hr 12 -2n,r 12 j (60) 

with n £E Z, while with 2j > non-contractible lines they are 

he 1 n+t2T2 — 1 — 2n,£iri+e 2 »'2 — l+2j ■ (61) 

The sign ei = ±1 indicate whether the leftmost non-contractible line is required to touch the 
left boundary (for ei = +1, referred to as the "blobbed" sector), or forbidden from doing so (for 
ei = —1, referred to as the "unblobbed" sector). The sign e 2 similarly describe the choice of 
blobbed/unblobbed sectors at the right boundary, and n £ N. The parameters r%, r 2 and T\ 2 are 
related to the weight of marked loops through the formulas (59). The case fa = (3 and fa 2 = fa 
corresponds to the one-boundary loop model (1BLM) [60], where the exponents reduce to [61] 

(62) 

5.1 One-boundary case 

We will start by investigating the one-boundary problem. Setting n — 0, the Hamiltonian we want 
to study is 

ra-l 2L+to-2 

H=-uJ2 p i,i+i- £ E i- (63) 

i— i— m 

When computing the spectrum of this operator, one extra representation theoretical consider- 
ation comes into play. Once the number of non-contractible lines (or "strings") at the boundary 
have been fixed, (63) still acts on them non-trivially by means of permutations. The Hamiltonian 
can be block-diagonalized with respect to this action by a change of basis that fixes the action of 
the symmetric group & m +j on the m + j strings according to its irreducible representations, the 
so-called Specht modules [62]. Since permutations occur between boundary lines and the leftmost 
bulk line (if any), we will actually deal with & m when j = or <5 m +i when j > 0. 

We can follow the procedure explained in [15] for a related model based on the Brauer algebra. 
Instead of the usual diagrams on which the Hamiltonian acts, we can consider a tensor product of an 
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unlabcllcd diagram, for which states with exchanged strings are equivalent, and a vector belonging 
to a basis of a certain Specht module. An operator will act on such a state by the multiplication of 
diagrams, and by the matrix action in the Specht module of the permutation keeping track of the 
reordering of strings. 

Then the Hamiltonian (63) can be put in a block form, each block indexed by the number of 
bulk strings and a Young diagram referring to a Specht module. Accordingly the Hilbert space 
will be decomposed in terms of these subspaces. Numerical diagonalization shows that the lowest 
eigenvalue of the Hamiltonian always lies in the symmetric representation of the symmetric group, 
where permutations act as the identity in the Specht module. So for the purpose of computing the 
leading critical exponents we can restrict to this sector. 

5.1.1 Exponents and relations with the one-boundary loop model 

Diagonalizing numerically the Hamiltonian at different values of /3 and to, we have been able to con- 
jecture the analytical form of the (m + 2j)-leg critical exponents h m -°(j). Using the parametrization 
of the 1BLM, if we define /3i as 

01 = ^TT ' (64) 



then the exponents are 



h m '°(j) — ft ri ( J ri), ri ( m )+2j ; (65) 



l±P ow f« ( (!+"») 



where r% (m) is given explicitly by 

n (m) = ^ arctan ( ^p^pSLl _ j , , , )( , , 

Our numerical findings give strong evidence that this result for the exponents does not depend on 
the coupling u as long as we take it positive (see figure 9). Note also that when (3=1 (whence 
Pi = fx (to) = 1 for all to), the exponents are trivial, equal to those of the free boundary case. 3 

In view of this result, which is formally identical to (62) for a particular choice of the parameter 
ri(m), it is natural to search for a relation between our model and the one-boundary Temperley- 
Lieb algebra. We recall that the latter is defined [60, 61] by endowing the Temperley-Lieb algebra 
with an extra blob operator b acting on the leftmost strand, which we will label to for making easier 
the connection with our model. By definition, b is idempotent and satisfies the relation 

E m bE m = (3iE m , (67) 



3 If we think of the case m = 1, one intuitive argument supporting this triviality goes as follows. Everytime a loop 
crosses the boundary line, we can think of this configuration as a loop with weight 1 and a straight line times the 
coupling — u: 



o 



For general m and u positive we have found the following relation = — u + AjL being the lowest eigenvalues 

of the Hamiltonian with m boundary lines, which implies the same scaling in 1/L of the two eigenvalues. 
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One-boundary (3 — 
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m = 1 , j = i 1 1 _ 

m = X, j = 1 h— -x — i 
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Figure 9: Independence of critical exponents h on the coupling u for j = 0, 1, 2 and m = 1, 2 (n = 0) 
in the model with j3 = 0. The values are computed by finite-size scaling corrections of equation 
(58) with inclusion of a term of order 1/L 2 , for systems of sizes L = 5 to L = 9. 

whose geometrical interpretation is that marked loops get a weight j3\ . 

For the particular case m — 1 this relationship is manifest. Namely, setting Pq = 26 — 1 and using 
b 2 = b we get (-Po) 2 = 1 as required. So for a particular choice of the parameter u the boundary 
interaction is proportional to b indeed. The corresponding weight of boundary loops follows from 

*i±fl- a ,68, 

where we have used (48). This shows that fix — -Mp as in (64). The exponents (65) then follow by 
invoking the result (62) for the one-boundary loop model [61]. 

For higher values m > 1 there is no such mapping to the one-boundary loop model. We can 
nevertheless derive the results (64)-(66) for an appropriate modification of the Hamiltonian (63). 
Our argument that (64)-(66) are correct also for the original Hamiltonian (63) is then based on 
universality: we present numerical evidence that the original and modified Hamiltonians have the 
same critical exponents. 

As a warmup, we first discuss our numerical observation that within the space of Hamiltonians 
of the form (63), the exponents (65) are independent of the boundary coupling u. This should be 
interpreted in terms of boundary renormalization group flow. The situation is akin to that of the 
Ising model with a boundary field, where the fixed boundary condition is a stable fixed point, and 
the free boundary condition an unstable one [63] . At a given non-zero value of the boundary field, 
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there is a crossover, and the system flows to the stable fixed point for a sufficiently large system. 
In our case, the boundary perturbation to the Lagrangian is u<j) rx r m \ ri r m ), and it drives the system 
from the m = to the m > boundary condition. Since the dimension of the field 4'r 1 (m),r 1 (m) 
is h ri f m \ ri f m ), by dimensional arguments the relevant length scale in this boundary flow is given 



by L 



1-hr 



We have explicitly checked that curves describing the critical exponents as a 



function of u, collapse after appropriate rescaling in the st(l|l) case, where we can access large sizes 
of the system, see figure 10. 
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Figure 10: h vs u in the s[(l|l) chain for sizes of the bulk length L = 45 to L — 75 (here in different 
colors) and m = l,n = 0. In figure (a) the variable u is not rescaled. In figure (b) we rescale u as 



L 



l-h r 



> r i( 1 )u, and we observe a perfect collapse of all curves. 



Extending these ideas of universality, we now construct a modified Hamiltonian H for the lattice 
model with edge states, which — unlike the original Hamiltonian H of equation (63) — is in precise 
correspondence with the one-boundary loop model. The key step is to construct an idempotent 
blob operator b satisfying (67). 

In the Hilbert space of our spin chain, the idempotent words of the boundary algebra are given 
by the Young symmetrizers of the symmetric group <3 m +i. In appendix B we show that indeed 
every symmetrizer satisfies a relation such as (67). The value of ft (64) we are interested in is 
obtained if we take the fully symmetric representation of the symmetric group (see figure 11 for a 
graphical representation of the relation between the symmetrizer and the blob operator). 

Define therefore the standard Young tableau 4 t\ — 1 I ' ' ' 



and the associated Young 



symmetrizer bt t = 



(m+l)! 



a. It then turns out that bt t satisfies: 



Emb^Em — 



m - 



-jE m b {i , 



(69) 



4 In what follows we will adopt the convention of labelling the boxes of a Young diagram of size m + l with the 
elements of the set {0, . . . , m} when related to the left boundary action of the permutation group 
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i 



Figure 11: Symmetrization (here represented by horizontal gray bars) over the to + 1 leftmost 
strands can be interpreted as blobbing the first bulk strand. 



where i\ = 1 • • • I to — 1 1 . To prove formula (69) we can use the following simple identity 

b tl = — — -b± (1 + (0, m) + (1, to) H h (m - 1, m)) , (70) 

TO + 1 1 

denoting the transposition of i and j, and the relations of the algebra. 
Relation (67) now follows with our value of /3\ (64) if we replace Ei by E% — Eib^ , an operation 
which does not affect relations of the algebra -4.2£,m,n- In other words, the following modified 
Hamiltonian 

j 2L+m-2 

H = -u- — V a- V Ei, (71) 

(to + 1! ^ ^— ' 

o"G©m + l i=m 

has precisely the form of a one-boundary loop model (the first term being —ub) and by (69) the 
weight of a loop marked by b is exactly (64). It follows from (62) that the 2j-leg exponents are 
indeed given by (65). 

Obviously the boundary terms in (63) and (71) are different. However, just as the spectrum 
of (63) was shown numerically to be independent of u, it is a natural hypothesis that the critical 
exponents of H and H are identical. In table l we present strong numerical support for this 
universality hypothesis. 

We remark that for having the universal behaviour described above, it is crucial that the signs 
in front of the permutation generators are all negative. In this case, adding to the Hamiltonian (63) 
other terms of the symmetric group than the generators, is a redundant perturbation when they 
have the same negative coefficients. We expect that taking different couplings Ui for each Pi.i+i in 
the boundary Hamiltonian, will in general drive the model to different universality classes according 
to the relative signs (see also appendix B). In particular, it appears an appealing possibility that 
each choice of Young symmetrizer in the above construction might lead to a different universality 
class, but we have not investigated this issue in sufficient detail. 

Finally we have numerical evidence that the universality class of H and H is more general, and 
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Tabic 1: Comparison of critical exponents for the model described by (63), h, and for the one 
described by (71), h, listed for different values of j3 and j, in a system with m — 2,n — 0. The 
values are obtained by fitting eigenvalues for bulk sizes L = 6 — > 10, using formula (58) with a 
term of order 1/L 2 , and although these numbers should be equal when L — > oo, already with these 
relatively small sizes we find good agreement between the values of the exponents in the model and 
the exact value predicted in the continuum limit /i2,oO)i given in the last column. 

that, in particular, the following Hamiltonian also has the same continuum limit: 

H = b {i Hb h , (72) 

corresponding to the case where the m strands on the boundary are fully symmetrized. From 
an algebraic point of view this version of the model is equivalent to replacing the m copies of 
the fundamental representation of the superalgebra in the spin chain formulation with a higher- 
dimensional one. 

5.2 Two-boundary case 

We now move to the two-boundary problem, which is obviously much richer. 
5.2.1 New sectors 

The quotient given by the action of the Hamiltonian on the super spin chain discussed previously in 
section 4.3 allows contraction of boundary lines, introducing a new quantum number in the problem: 
k, the number of couples of left-right boundary lines not contracted. So we will restrict to (J, k) 
sectors with j and k fixed and the action of TL generators on a pair of strings (non contractible 
lines) set to zero. See figure 12 for examples of states in different sectors. 
A sector of given k appears with multiplicity 
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Figure 12: Different sectors in a system with m = n = 2 and L = 3. From left to right, (j = 0, fc = 
0), (i = 0,fc = l), (j = 2,fc = 2). 

where recall that m = min(m, n). The possible values of k for j = 0, 1 are k = 0, 1, . . . , m. Indeed 
when j = or j = 1 we can always go from a configuration with k > contracted boundary lines 
to another with fc — 1 contracted boundary lines by acting with elements of the algebra. In figure 
13, it is shown how we can go from the sector fc = mtofc = m— 1, and it is clear that previous 
exchange of boundary lines allows contracting up to m lines. For j > 1 meanwhile, there is only 
the sector k = m. In this case the boundary lines cannot come on adjacent sites since lines in the 
bulk cannot cross, as illustrated for j = 2. 

Note that with this definition, the sectors (j — 0, k) and (j = 1, k — 1) are actually the same, 
and that a state can be regarded as belonging to one or the other sector (see figure 14). 

Therefore it remains to study the sectors (j = 0, k = 0, 1, . . . , m) and (j, k = m) when j > 0. In 
each such sector the total (bulk plus boundary) number of non-contractible lines is different and 
correspondingly the value of the exponents will depend on j + k. Calling £ — \m — n\, we will refer 
to the lowest exponent in a sector (j, k) as the (2j + 2k+£)-leg exponent, h m ' n (j + k). Note anyway 
that if we would like to look at the lowest exponent when we fix only j, this will be present in the 
sector with the minimum value of k possible. 

5.2.2 Relations with two-boundary loop models and solution 

In the analysis of the one-boundary problem we have found that our model is in the same universality 
class as a 1BLM with weight of marked loops given by (64). This weight can be obtained by taking 
as blob operator the Young symmetrizer of the fully symmetric representation of the symmetric 
group. This identification now turns out to be very useful for solving the two-boundary problem 
which we recall is described by 

m-l 2L+m-2 2L+m+n-2 

H = -uJ2 Pi,i+i ~ Ei ~ v J2 P M+i- (74) 

i—0 i—m i—2L J r ni—l 
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Figure 13: (a) Starting from the reference (reduced) state on top left at j = with fc = m, we pass 
to the state with k = m — 1 on the right upon acting with the word juxtaposed above, (b) Similarly, 
when j = 1 multiplication by an element of the algebra allows passing from fc = mtofc = m-l. (c) 
When j = 2, we have to exchange bulk strings in order that the two boundary strings (the dashed 
lines) can be permuted with the bulk strands, and come on adjacent sites. This is forbidden by the 
interactions in our model. The situation is similar with higher values of j. 
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Figure 14: From the reference state of the sector (j = 0, k) on top, we can reach the state of the 
sector (j = 1, k — 1) on bottom, by the two steps depicted. First, for each boundary we pass one of 
k lines into the bulk by applying permutations (middle), and then we order the strings using TL 
generators (we have assumed m > n). Going from bottom to top works in the same manner, and 
this procedure can be done for every states in the two sectors, which are then equivalent. 
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Indeed, by naive analogy, we expect the following Hamiltonian to be in the same universality class 
as our starting problem (74): 

| 2L+m— 2 

H = -u- ' " 

( //; - I i j 

-Ub tl - ; c; - r <),, 



j 2L+m-2 ^ 

(mTT)! S ff E ^- U (^TT)! E CT 

L+m-2 

E Si-w6t 



(75) 
(76) 



where we have introduced the symmetric group of the m + 1 leftmost strands £> lctt = @/o,...,to}i 
the symmetric group of the n + 1 rightmost strands £> rlght = @{2i+m-i,...,2i+m 
the standard Young tableaux t\ = 1 I ' ' ' M and i2 = I 2L + 777 — 1 | ■ ■ ■ | 2L - 



_n, as well as 
n + ti - 1 I As 



before 6 t stands for the Young symmetrizer related to the tableau t. 

We have checked numerically the hypothesis that adding these additional permutations to our 
former Hamiltonian H is a redundant perturbation also in this two-boundary case. See table 2 for 
a sample of evidence in the case m = n = 2. As in the one-boundary case, we claim that the results 
are unchanged as long as u and v are taken positive, as we have verified this numerically (see figure 
15). 
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0.3346 ± 0.0003 


0.3346 ±0.0003 
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2.118 ±0.008 


2.118 ±0.008 
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Table 2: Comparison of critical exponents for the model described by (74), h, and for the one 
described by (76), h, listed for different values of j3, j and k, in a system with m — 2,n = 2. The 
values are obtained by fitting eigenvalues for bulk sizes L = 4 — > 8, and although these numbers 
should be equal when L — > oo, already with these relatively small sizes we have good agreement 
between these two values and with the exact value predicted, given in the last column. 

The advantage of studying H is that it belongs to the two-boundary Temperley-Lieb algebra, 
the two-boundary analogue of the blob algebra, and for that problem we know how to compute the 
exponents [59]. In the language of the two-boundary loop model (2BLM), loops marked on the left 
boundary have weight (64) and those marked on the right boundary 

i + /3 



02 



1 



(77) 



5 CRITICAL EXPONENTS 



31 



Two-boundaries /3 = 1 



0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 


-0.05 



1 

• • 


i 

• -m- 


i 

-* — * 


i i 


i 

--•—»—«- 


1 

-*—-«- 


rn 


= 71 = 




= 0,fc 


= i — 




m 


= n = 


i,i 


= 0,fc 


= 1 - 


--x 1 


m 


= 71 = 


i,j 


= l,fc 


= 1 - 


: 


m = 


2,n = 




= 0,k 


= 


■-B -. 


m = 


2,n = 


i,i 


= 0,k 


= 1 h- 


■ 1 


m = 


2,n = 




= l,fc 


= 1 h. ~ 





SF--K— *e— x— x— ^- 



■x — — k — x — -x- — x--x — — k — a — n — x--^ — 5! 



1 



4 5 6 

u = v 



10 



Figure 15: Independence of the first few critical exponents on u = v for m = n = 1 and m = 2, n = 1 
when (8 = 1. The values are computed by finite size scaling fitting equation (58) with a term of 
order 1/L 2 , for systems of sizes L = 4 to L = 8. 



All that remains to be done is to compute /3i2, the value of loops marked by both blob operators, 
as expressed in the following algebraic relation (we can consider for simplicity the situation with 
L = 1; in the general case replace E m by E defined in (50)): 

E m bib2E m = /3i2-E m ) (78) 

where b\ and 62 act respectively on the first and second strand labelled m and m + 1. 

So wc fix L = 1, and consider the j = sector only. Define £1 = 1 • • • \m — 1 1 ( £ 2 = 
m + 2| ■■■ I m + rt + If we replace respectively 61 and 62 in equation (78) with our operators 
bt 1 and bt % , and -Eo with £" m = 25 m £>£ &£ , we would like that the following equation holds for a 
certain value of P12: 

E m b tl b t2 E m = (3 12 E m . (79) 

This equation can however only be satisfied and meaningful in our model, if we consider it for a 
fixed value of k, as we will explain below. The sector (j = 0, k) of the Hamiltonian (76) will then be 
described for different k's by different quotients fixing /3i%, which we will now call /3^ 2 to indicate the 
fc-dependence. Taking the quotient as described in section 4.3, we have already implicitly chosen 
particular values of /3f 2 i which we are going to compute. We remark that the representation of the 
two-boundary Temperley-Lieb algebra given by (76) is not faithful in general, and the structure of 
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the sectors cannot be directly given from this mapping. However the leading eigenvalues will be 
present and we explain now how to compute them in the different sectors. 

Before showing the general procedure, we will give an example in the case m = n = 2, where 
we can understand how the mechanism works. First it is convenient to express b tl and b t2 using 
the following identities: 



\b i% (1 + (0, 2) + (1,2)) - \b h (1 + Pi(l + Pit)) 



6 t2 =^(l + (3,4) + (3,5))6 f2 



1 



(! + (! + P 4 )P 3 )& t - 



(80) 
(81) 



Then, using the commutation relations of the algebra, we have: 

E 2 b tl b t2 E 2 = ^E 2 b k (l + Pi(l + P ))(l + (1 + P 4 )P 3 )b t2 E 2 
= ~ (08 + 4)P 2 + 4& fi & f3 BaPiiWfX) ■ 



(82) 



The first piece of equation (82) is already in the desired form, while the second piece is not propor- 
tional to E%. As anticipated, we can only achieve this by fixing a sector k, and computing the action 
of the operator W :— bf b^E^PxPzE^bi b^ on the states belonging to the sector. The graphical 
representation of W and E 2 is depicted in figure 16. 



W = 




E 2 = 



Figure 16: The words W and E 2 in the diagram representation. As before, the gray bars indicate 
symmctrization over the lines and W is obtained using the quotient adopted in section 4.4. 



Let us first consider the case k = 2, when no contractions of boundary lines is allowed. It is 
clear from figure 16 that W acts as zero on every state in this sector, since it diminishes the number 
of strings from 2 to 1. So, we have found that every state of the sector k — 2 is in the kernel of the 
operator 

R := (E 2 b tl b t2 - ft 



Ei 



(83) 



with 



/3i 2 2 = 



13 + 4 
9 



(84) 



We now move to k = 1. In this sector we have 18 states. They are drawn in figure 17, where 
we have grouped those related by the action of b^ b^ = (1 + Pq + Pi + PqP\)/A. Acting with E 2 on 
an arbitrary state will always produce the sum of the four states contained in the box in figure 17 
times a constant which can eventually be zero, as when E 2 acts on two strings. The role of E 2 is 
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Figure 17: The 18 reduced states of the sector k = 1, for a system with m = n = 2, L = 1. 
related by the action of are grouped together. 
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indeed to project onto the state in which the boundary strands are fully symmetrized and the bulk 
ones are contracted. This remark is of crucial importance since it tells us that the action of R is 
zero on every state, with the same value of /3f 2 then characterizing this sector. We can evaluate R 
on a reference state, say the one on the top left in figure 17 to obtain the value of /3f 2 - R on this 
state is zero if: 

2/3- 



#2 = 



6 



9 



(85) 



The last case we have to consider is k — 0, when no strings are left. In this sector we have 6 
states depicted in figure 18, where again we have grouped those connected by the action of b^b^. 
As in the previous case, we realize that Ei acting on every such states produces the sum of the two 





Figure 18: The 6 reduced states of the sector k = 0, for a system with m = n = 2, L — 1. States 
on top and on bottom grouped together are related by the action of b$ b^ . 

states contained in the box in figure 18 times a constant. Then, evaluating R on a given state gives 

13 + 2 



(86) 



Having understood this simple example, we can consider the general case with m and n arbitrary. 
We first express b tl and bt 2 in terms of b$ and 6 t ~ 2 : 



bu - 



— bf (1 + (0,m) + h (m — l,m)) 

- 1 

b t2 = ~~T^ (1 + (m + 1, m + 2) H h (m + 1, m + n + l))b t 



ir( i +(En p ™ + ^) 

\ \i=0j=0 / 



Pm+1 bi 



(87) 



5 CRITICAL EXPONENTS 



35 



Then wc arrive at the generalization of formula (82): 



E m b tl b t2 E m — 



(to + l)(n + 1) 



1 



( 



(/3 + to + n)£ m + 



(89) 



(to + l)(n + 1) 



\ j=0 j=0 

1 



/ n — 1 i— 1 




) , (91) 



(90) 



where the last equality holds thanks to the fact that £>£ can be written as itself times the Young 



symmetrizer related to the tableau obtained by removing the box numbered m — 1 in t%, and the 
analogous one for 6 i2 . 

Now, we can generalize our reasoning done for the case m — 2, n — 2. If we redefine the operator 



then every state of a sector with given k will be in the kernel of R with the same value of /3f 2 i a 
consequence of the projection performed by E m . The values of /3f 2 m the general case to, ti are 
then found to be: 



Note that (3\ 2 1S symmetric in exchanging to, n, as it should in our problem. The computation of 
the weights /3f 2 m the model described by H, equation (76), concludes the analysis of the sector 



As explained in section 5.2.1, when we fix j > 1, we have only to consider k = m, because if 
j = 1, the other values of k are related to corresponding sectors in the j — case, and when j > 1, 
no contraction of strings is allowed. So for j > 1, the mapping onto the 2BLM can be directly done 
without any further issue. 

The (2j + 2k + £)-\eg exponents h m , n (j + k) in the model (76) are finally listed in table 3 where 
we have used the parametrization of /3i , /32 , /3i2 m terms of r\,r2,T\2 introduced in (59). Recall 
that j3\ and f}% are simple functions of j3 and of respectively m and n, equations (64) and (77), and 
that m is the minimum between to and n. When the number of legs is between t and to + n, we 
are in the j — sector of a 2BLM, and for k = 0, . . . , m, the exponents h r k^ r k^ are computed from 
our expression of /3f 2 , equation (93). When we have instead m + n + 2j legs, with j > 0, we are in 
the sector j of a 2BLM which is completely described by the two parameters @x and /?2, equations 
(64) and (77), and the exponents then follow from fusion of the boundary conformal weights of the 
one-boundary case, the leading exponents in a sector being ft. r i+r 2 -i.ri+r 2 -i+2j- 




(92) 




(93) 
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Table 3: Critical exponents for the two-boundary problem. 

Invoking the usual universality hypothesis — namely that the critical behaviour is unchanged 
upon adding permutations other than the generators to the Hamiltonian H — the above expressions 
completely determine the critical exponents in our original model (74). 

Lastly, we have verified that as in the one-boundary case, the version of the model with higher 
spins on the boundary, H = b$ 2 Hb^ b$ , is in the same universality class as well. 

6 Back to the spin chains spectrum 

In section 4.3 we have suggested that the representation of the quotient algebra Q,2L,m,n{M) is 
faithful if the number of states per site in the spin chain is greater or equal to three. Therefore, 
every irreducible representation of the diagram algebra studied previously should be present in the 
sl(N + M\N) spin chains, except when N = 1, M = 0. Whenever the representation is faithful, 
the eigenvalues of the loop and of the super spin representations are the same, the difference in the 
Hilbert spaces manifesting itself only in the degeneracies. We can thus use the loop model — which 
is technically more convenient — to solve the full spin chain problem. This requires understanding 
all exponents, and calculating their degeneracies. While such a program can easily be carried out 
for the pure Temperley-Lieb problem [17], in the present case, a lot is still lacking. 

One first incomplete aspect of our study is that we have not determined all the geometrical 
exponents: there are definitely plenty of subleading eigenvalues in a given sector. When a direct 
mapping onto a Boundary Loop Model is available, we can reinterpret the (J, k) sectors of the edge 
state model in terms of the blobbed/unblobbed sectors of the BLM and we hence know exactly 
which conformal weights one finds. Such an exact mapping is possible in the one-boundary system 
with m = 1 and m — 0, for which the conformal weights of primary operators in the sector j = 
are h ri , ri , and for j > are h ri ^ ri+ 2j (blobbed sector) and /i r i,n-2j (unblobbed sector). It is 
also possible in the two-boundary case with m = n = 1 if we fix k = 1, where for j = we find 
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primaries with weight fo r i 2 _2n,r 1 and n € Z, while for j > we find /i eiri + e2r2 — i-2«,eiri+e 2 r 2 -i+2j 
with t\ t 2 — ±1 and neff. But already in the sector j = k = when m = n = 1, our algebra gives 
an unfaithful representation of the two-boundary Temperley-Lieb algebra, since the mapping to a 
2BLM involves a quotient version of the latter, and we do not find all the subleading eigenvalues 
h r o 2n.r1 2 1 n £ Z. In general we do not have full control of which subleading eigenvalues appear. 
However the value of possible eigenvalues we expect should be of the form given by the 2BLM 
conformal weights with some possible omissions and blobbed/unblobbed sectors mixing, as we have 
indeed observed in the systems studied numerically. Clearly a deeper algebraic understanding of 
the problem is needed here to characterize the full spectrum. 

The other incomplete aspect is that we lack a clear understanding of the degeneracies. One 
way to understand these in general is to think of the generating function of eigenvalues (for the 
super spin chain), which could in principle be computed as the modified partition function of a 
loop model where non-contractible loops are weighted differently [17]. In algebraic terms, this 
generating function would be written as a sum over irreducible representations of the algebra to 
which the Hamiltonian belongs. The multiplicities of summands would then be the dimensions of 
irreducibles of the commutant algebra in the Hilbert space at hand; the problem is that very little 
is known in general about this commutant. 

Another way to attack the problem of degeneracies is to study the combinatorial aspect of the 
loop representation of the model, along the lines of [64]. One would then need to develop a theory 
of combined projectors of the symmetric group (on the boundary) and of the Jones- Wenzl type (in 
the bulk) in which the degeneracies would be identified with the Markov traces of these projectors. 
One way or the other, it is clear that further study is needed to finish the program. 

For some physical applications [32], only the leading exponents will matter, and the exact 
knowledge of degeneracies and subleading corrections to scaling does not play a paramount role. 
However, if these two aspects were completely elucidated, one could construct exact continuum 
limit partition functions and compute various kinds of exact crossing probabilities, presumably with 
applications to transport properties in the spin quantum Hall effect. We leave these developments 
for future work. 

Let us finally recall that the s[(l|l) super spin chain — alias symplectic fermions; see appendix 
A — seems to be the only case where the spin chain representation is not faithful. In particular this 
manifests itself in the fact that for m — n, the spectrum in this chain is equal to that obtained 
for free (m = n = 0) boundary conditions, when the Hamiltonian belongs to the Temperley-Lieb 
algebra and the critical exponents are those of the conformal field theory of symplectic fermions. 
Instead, as explicitly showed in table 4, in the case m — n ^ 0, the sI(A|A) spin chain when N > 1 
possesses an infinity of new exponents which are not contained in symplectic fermion theory. Note 
that these exponents in general are not rational — a rather unusual feature in conformal field theory. 
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Critical Exponents (3 
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Table 4: In the first column, watermelon exponents computed in the geometrical model when 
m = n = 1 and /3 = 0. The number of legs is the coordination of the vertex. In the second and 
third column we note the presence or not of these exponents in the faithful (N > 1) and unfaithful 
(N = 1) representations of the algebra for the sl(N\N) spin chains. 

7 Conclusion 

In conclusion, we have seen that edge states in super spin chains describe an infinity of conformal 
invariant boundary conditions for conformal superprojective sigma models. The technology involved 
to characterize these boundary conditions is definitely involved, and the exponents very far from 
trivial. How these exponents can be organized in a complete description of the boundary and bulk 
logarithmic CFT describing the low energy physics of the sigma models remains an open question. 

Meanwhile, it is natural to expect that the boundary conditions we have uncovered, together 
with their edge states description in the quantum spin chains, will play an important role in the 
description of transitions between plateaus in the spin quantum Hall effect, as well as in the ordinary 
quantum Hall effect (though in this case a different sigma model would have to be considered) [29] . 
We hope to get back to this question soon. 
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A The spin chain 

Here we will work out in details the solutions of the model in the simplest case of the flt(l|l) super 
spin chain, where the Hamiltonian reduces to free fermions and explicit computations can be done. 
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A.l The lie superalgebra gl(l|l) 

First we will briefly recall some facts about the Lie superalgebra g[(l|l) and its representation 
theory [65]. 

The g[(l|l) Lie superalgebra is defined by two bosonic generators E,N and two fermionic ones 
F~ , F + , and its superdimension is zero. E is central, and the other defining relations are 

[N,F ± ] = ±F ± , {F-,F+} = E. (94) 

In particular N counts the number of fermions. 

When dealing with the irreducible representations of superalgebras, we have to distinguish be- 
tween typical or atypical representations. For |1) typical representations are the two-dimensional 
representations (e, n), e 7^ 0, n € R, which read in matrix notation: 

-(V I) F = (I I) r>-(l I) o», 

The state annihilated by F~ is a boson. Atypical representations are the one-dimensional irreducible 
representations (n), which have F + = F~ = E = 0, N = n. 

Other finite-dimensional representations that occur are indecomposable, not fully reducible. 
They are composites of atypical constituents linked by the action of F ± . We denote by V n the 
indecomposable four-dimensional representation playing the role of the projective cover of (n), that 
is it can never appear as a subrepresentation of a larger indecomposable and contains (n) . 

The spin chain (37) is built with the fundamental V = (1, 1/2) and the dual V* = (—1, 1/2). 
The (graded) tensor product of representations is given by the following formulas: 



(ei,ni> <X> (e 2 ,n 2 > = 



£ni+n 2 -i ifei + e 2 = 0, 

(ei + e 2 , ni + n 2 ) © (ex + e 2 , m + n 2 — 1) otherwise. 

(e, n) (g) V m = (e, n + m + 1) ® 2 • (e, n + m) ffi (e, n + m - 1) . (96) 



Then we have for l± ^ Z 2 : 

V® h 9 (V*)«« - © © © f 1 7 ( l 2 7 X ) (h + h, i + J + 2 - ^ + p) . (97) 

This is the decomposition of the Hilbert space (37) when the number of extra representations on 
the two boundaries of the spin chain is different. It involves jjl(l|l) 2-dimensional representations, 
since changing the order of representations in the spin chain leads to isomorphic spaces. When 
instead l± = I2 = I 

11 '21-2 



{V ® V*)® 1 ~ © [._ * )Vi, (98) 
i=l-i ^ ' 

this corresponding to the decomposition of the Hilbert space in terms of 4-dimensional representa- 
tions when we have the same number of representations at the two boundaries of the chain. 
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A. 2 Free fermions formulation 

Denote by fi,fj the fermionic operators of 1( 1 1 1 ) acting on V at site i, and by fi,f} those acting 
on V* . They satisfy the anti-commutation relations {/•£,/!} = Sij, {fi, fj} = —Sij, where the 
presence of the additional minus sign implies that the Hilbert space has indefinite inner product. 
To deal with ordinary fermions Ci,c\ that satisfy {ci,cj} = Sij, we define Cj — fi and cj = fj if we 
have V on site i; if we have V* on site i define instead Ci — fi and c- = — fj. 
Pi and Ei defined in section 4.3 can be written in terms of these fermions as: 

Pjj+i = (1 - {c]c 3 - c) +1 Cj - ctcj+i + ct +1 Cj+i)), i = 0, . . . , m - 1 (99) 

Ej = (-iy +m+1 (c] Cj - c] +1 c 3 + c)c i+1 - c] +1 c J+1 ), j = m, . . . ,2L + m — 2 (100) 

Pj,j+i = -(1 - (c]cj - c] +1 Cj - cjcj+i + c] +1 Cj+i)), j = 2i + m - 1, . . . , 2L + m + n - 2 . 

(101) 

One can check by hand that the relations (46)-(49) are satisfied. The Hamiltonian (40) in this 
representation becomes quadratic in the fermionic operators: 

m 1 

H = u ( c ] c j - c )+i°j - c ] c i+i + c ]+i c ]+i)+ 



3=0 

2L+m-2 

E ( 

3=m 



-l) j+m (c] Cj - ct +1 cj + c]c,+i - ct +1 c i+1 )+ (102) 



2L+m+n-2 

-u ^ (ctcj-ct +1 Cj-ctc j+ i + ct +1 Cj + i)-mu + ra;. 

j'=2L+m-l 

A. 2.1 Diagonalization procedure 

Rewrite (102) as 

H = cjAijCj — mu + nv , (103) 

i,3 

where A is an appropriate (2L + m + n) x (2L + m + n) real non-symmetric tridiagonal matrix. 

The problem of diagonalizing a Hamiltonian of the form (103) is a well-known problem, solved 
by Lieb, Mattis and Schultz in [66]. We introduce by a canonical transformation the fermions Tjk,m., 
in terms of which H is diagonal: 

H = J2^(vUk-iy (104) 

However the equations involved in the diagonalization procedure even in the simplest case 
m = l,7i = are rather complicated, and instead we will proceed diagonalizing the matrix nu- 
merically This task we can accomplish more efficiently by using a closed formula for the charac- 
teristic polynomial of A for generic values of m and n. Indeed we have found that P L '"(it, v, x) — 
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T)et(A — xl2L+m+n) has the following expression: 
P™' n (u,v,x) = x(-l) L+n u™- 1 v n - 1 ! { U L+1 

U m +i (l - 7^) (pi(u,v,x)U n (l + -^j +p 2 (u,v,x)U n+1 (1 + -^j ) 
+ U m (l - (p 3 (u,v,x)U n (l + -^-j +pi(-v,u,-x)U n+ i (l + ^- 



u L \i-- 



U m +i (l - ^) (qi{u,v,x)U n (l + +g 2 (u,w,x)[/ n+ i (i+ |^)) + 
+ U m (l - (q 3 (u,v,x)U n (l + -^j +q 1 (-v,u,-x)U n+ i (l + ^ 



(105) 



where Ui(z) are the Chebyshev polynomials of the second kind, 

, sin ( (i + 1) cos _1 (z)) 
Ui{z) = U K ,} , (106) 

V 1 — z z 

and the coefficients p's and q's are simple polynomials depending only on u, v and x: 

Pi(u, v,x)=u (x 4 + 2x 3 v - x 2 (v + 2) - Axv + v) (107) 

P2{u, v, x) = —x (x 2 — 2) uv (108) 

p 3 (u,v,x) = x (x 4 + 2x 3 (v - u) - x 2 (4:uv + u + v + 2) + 4x(u — v) + 9uv + u + v) (109) 

qi(u,v,x) =-u(v-x (x 2 - 3) (x 3 + 2x 2 v - x(v + 1) - 2v)) (110) 

q%{u, v, x) — —x (x 4 — Ax 2 + 3) uv (HI) 
q 3 (u, u, x) = — x(t> (a; (a; (x (4ux — 2x 2 + x + 8) — 17u — 3) — 6) + 14u) 

+ x (x 2 - 3) (x (-2ux + u + x 2 - l) + 2u) + u + v) (112) 

Formula (105) seems to not help in the analytic determination of the eigenvalues apart from 
the trivial case m — 0,n = 0, but it is remarkable that the dependence on L, m, n is organized 
in terms of Chebyshev polynomials, and further it allows efficient numerical diagonalization of the 
Hamiltonian. The presence of Chebyshev polynomials does not come as a surprise since it is an old 
result that these polynomials are related to the discrete Laplacian [67], and our formula is in this 
sense a generalization of that result. 

We remark that the matrix A introduced above is actually not diagonalizable when m = n, and 
it has Jordan cells of rank 2. This is linked to the indecomposable nature of the conformal field 
theory describing the chain [19, 25]. 
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A. 3 Spectrum of the Hamiltonian and continuum limit 

We would like now to determine the spectrum of the Hamiltonian (104). We define the vacuum |0) 
as the state with no particles: % |0) = for every k, whose energy is — | Afc. Define I" and I + 
as the set of values of fc's such that the real part of the corresponding eigenvalues are respectively 
negative and positive. The ground state |f2) is obtained by filling the single particle states with 
negative energy up to the Fermi surface, = HfeeJ- 1^)' anc ^ ^ s ener gy ls 



^=Me a *-e a ■ ( n3 ) 

\kei- kei+ I 

We can calculate the conformal weights hi from finite-size effects using formula (58). Ultimately 
with the information gained from finite-size scaling, we can compute for given boundary conditions 
m, n: 

Z mn = lim TYe -P(H-E (L)L) = TTq L -c/24 = y q h t -c/2^ (m) 

i 

where Tr stands for an ordinary trace on the vector space, not a supertrace, and q is as usual the 
modular parameter. This corresponds to a modified partition function, with antiperiodic boundary 
conditions in the time direction [17]. 

A. 3.1 Free boundary conditions 

We start by analyzing the simplest and already known case of the Temperley-Lieb chain with open 
boundary conditions (m = n = 0). In this case indeed formula (105) reduces to 

P° L '°(x) = , (115) 



and its roots are given by x — with multiplicity two, and the roots of Ul-i ( 1 — \ 



x = 2 ^ n [jlj > k = ±l,...,±(L-l). (116) 

The eigenvectors corresponding to these eigenvalues can also be easily computed. 

With these values of determined, we can use the Euler-Maclaurin formula and equation (58) 
to verify that the central charge of the model is c = —2 (using that v s — 2 for p = 1): 



Eo 1 vL . / k \ l-cot(^) 

— = > 2sm 7r — 1 — 

2L 2L ^ \ 2L 

fe=i v 

2 1 7T 

~ ~^ ~ 2L + 6(2L) 2 



(117) 
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To compute the partition function in the free boundary case, we note first that ho = 0. Then 
since £/ 7I "Si=i ^fc; ~ Si=i f° r l ar g e L, p L, the exponents are all integers and their com- 
putation is rather simple. Indeed it boils down to compute all the ways we can build the integer 
5Z?=i hi using at most two equal ki's (one associated to the action of ct., the other of c_fe 4 ). If we 
write 

oo 
i=0 

and use the notation A = (A" 1 , . . . , XZ k ) h i for the partition of size i, a, is given by 

fc 



4 ^ H 2 ^' 1 - ( 119 ) 



AH i=l 
n 3 =l,2 

The factor 4 comes from having two zero Afc's and parts occurring only once have a weight two 
because they can come from the action of ct. or c_fc 4 . As customary when dealing with the gener- 
ating function of partitions, we are not able to give an expression in closed form for aj but we can 
do this for the generating function: 

Zo.o = II i 1 + V 1 + Q 2h ) = 2^ = det(-V AtN ) , (120) 

h=0 ^^ T ' 

where T>a,n is the Laplacian on the strip with Neumann boundary conditions in the space direction 
and antiperiodic boundary conditions in the time direction. 

Finally we note that this result is in agreement with the expression of the partition function as 
sums of Virasoro characters (recall however that the modules involved are not simple but indecom- 
posable for this theory [19]): 

Z 0) o = £(2j + l) q - (121) 

= 2^ + 2y^ f| . (122) 

A. 3. 2 Generic boundary conditions 

We now turn to the discussion of generic values of m,n. In this case we do not have an analytic 
expression for the eigenvalues of the Hamiltonian and so we proceed numerically using the results 
presented in A. 2.1, and we extract the finite-size scaling by fitting energies for different sizes of the 
system. We have studied chains up to size 2L = 24. 
Define: 

2 / \m-n\ \ , 
r = 1 - - arctan 1 l — , 123) 
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where the absolute value comes from the symmetry upon exchanging m and n. We find that the 
partition function of the g[(l|l) chain is given by replacing the exponents /ii,i+fc by h r , r +k in (122): 



Now we will discuss the relation between this spin chain and the geometrical model at j3 = 0. 
We have pointed out in section 4.4 and 6 that the g[(l|l) representation is not faithful, so not all 
the eigenvalues of the Hamiltonian in the geometrical formulation are present here. We now wish 
to make this observation more quantitative. 

In the g[(l|l) spin chain the conformal weights appearing in (124) are h Tjr ±2j,j = 0, 1, . . . . If we 
compute explicitly (3\ 2 from our expression of exponents in the geometrical model (table 3) when 
P = 0, we find that the function r defined in (123) coincides with r^, so we interpret h r>r as the 
Meg exponent. Further r\ 2 = 2 — r and then ft, r>r +2 is the (£ + 2)-leg exponent. When we look at 
other exponents the situation is different for the one-boundary and the two-boundary cases. 

For the one-boundary case n = 0, in the spin chain we always find the eigenvalues of the 
Hamiltonian corresponding to the fully symmetric (and also to the alternating) representation of 
the symmetric group acting onto the strings, where the lowest eigenvalues lie. So, h r ^ r+ 2j are the 
(m + 2j)-legs critical exponents we found in the loop representation. 

When we look at the two-boundary problem, the lowest eigenvalues in the geometrical model 
for (j = 0, k = 0) and (j = 0, k = 1) are present in the spin chain, as already noted. However, 
for the other sectors this is no more the case, and the conformal weights /i r ,r±2j)i > 1) appear as 
subleading eigenvalues in the sectors of the geometrical model. 

We remark that for m — n the spectrum in the scaling limit of the g[(l|l) spin chain is the same 
as that in the free boundary case, and the multiplicity of the exponents h r , r ±2j is 4 while for m ^ n 
is 2. This can be interpreted from the analysis of the tensor product of representations of fll(l|l) 
given in (A.l). These numbers correspond to the multiplicities of modules of the commutant of 
gl(l|l), the Walled Brauer algebra, and are present also in our case. 

As a last remark, we interpret this result in terms of the continuum calculation presented in 
section 2.3 for the CP ' 1 sigma model, which the spin chain introduced in this section should 
discretize. Indeed we note that when [i = formula (19) can be rewritten as equation (124), where 
r = 1 - 2A. 

B Young symmetrizers as blob operators 

Irreducible representations of the symmetric group <3 m +i, the Specht modules, are obtained by a 
corresponding Young symmetrizer, indexed by a standard Young tableau t of size m + 1, where in 
each box of the diagram there are numbers from the set {0, . . . ,m}. Young symmetrizers b± have 




(124) 
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the property of being idempotent if properly normalized, b\ = bt- More than that, they satisfy also 
the other relation of the blob algebra: 

E m b t E m =p t 1 E m b i , (125) 

if we replace E m with E m = E rn bf. t is a standard Young tableau of size m, obtained from t by 
removing the box labelled m. So, one way to realize the blob algebra is taking a Young symmetrizer 
acting on the m boundary lines and the leftmost bulk strand. The weight of blobbed loops in this 
construction will be indexed by the Young tableau and is given by: 

FT r«A hook(cc) 
ILeY(t) book(a;) 

Y(t) is the Young diagrams to which the tableau t is associated, hook(:r) is the hook length of the 
box cc, Lt(m) is the number of boxes to the left of the box numbered m, and At(m) is the number 
of boxes above the box numbered m. 

One can verify formulas (125)-(126) by computing the first cases, doable by hand. For example, 
for the case of the Young symmetrizer of the symmetric representation, our formula is consistent 
with the result of equation (64). If we take the alternating representation instead, we have /3\ — 
((3 — m)/(m + 1). When one looks at 1BLM realizations of this algebra, however the interpretation 
of the critical behaviour of such a model can be difficult for negative f3\ , since phenomena like level 
crossings appear, as has been discussed in the context of the boundary chromatic polynomial [68]. 

We defer to future work further comments about the generalizations of this construction to the 
case of two-boundary Temperley-Lieb, which is definitely a more involved task. 
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